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THE EVOLUTION OF ELECTRON DENSITY AND TEMPERATURE 
DISTRIBUTIONS IN THE TOPSIDE IONOSPHERE DURING 
MAGNETIC STORMS* 

K. D. Coie ** 

J. A. Findlay 

ABSTRACT 

During geomagnetic storms the latitudinal distributions of electron density 
and temperature in the mid-latitude topside ionosphere are observed to change 
in a manner than can be related to the evolution of ring current particle popula- 
tions as predicted theoretically earlier (Cole 1964; 1965). The region of auroral 
precipitation is characterized by correlated increases in electron temperature 
and density. Equatorwards of this region there is a broad belt of elevated elec- 
tron temperatures and depressed electron densities which is usually much broader 
than any SAR-arc distinguishable from the ground, but which is nevertheless the 
same basic physical phenomenon. The changes of position of this belt can be 
related to prior bursts of geomagnetic activity and injection of ring current 
particles into the magnetosphere. 
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THE EVOLUTION OF ELECTRON DENSITY AND TEMPERATURE 
DISTRIBUTIONS IN THE TOPSIDE IONOSPHERE DURING 
MAGNETIC STORMS 

INTRODUCTION - , ; 

During the onset of the main phase of a geomagnetic storm electric fields 
associated with magnetic fluctuations pump plasma from the outer, to the 
inner magnetosphere. This process energizes. plasma to form a ring current 
{Cole 1964). Frank (1967) has observed fluxes of energetic protons (up to 40kev) 
in the magnetosphere sufficient to aceoimt for naost of the ring current during two 
geomagnetic storms in 1966. More recently energetic 0+, ions have been ob- 
served in the upper atmosphere at altitudes of about 800 km during geomagnetic 
storms over an L-value range of 2. 4 to 9 (Shelley et al. , 1972). Their fluxes 
are such that it appears that trapped 0+ ions could contribute significantly to the 
energy of the ring current. Cole (1970 a,b) has shown a correlation consistent 
with earlier theory (Cole 1964, 1965) between the latitudinal distribution, of geo- 
magnetic fluctuations and the positions of the ring current and the plasmapause 
in the cold plasma. A SAR-arc observed from- the ground during the storm 
of September 29, 1967 was related to the innermost L-shells of the relatively 
disturbed part of the magnetosphere. Figure 1 shows schematically a compari- 
son of groxmd based measurements of magnetic fluctuations during the onset 
of the main phase of a magnetic storm, proton fluxes, and plasmapause 
during the recovery of the main phase as deduced earlier (Cole 1970 a,b). 
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The association between the magnetic fluctuations and the location of ring cur- 
rent particles during the storm is employed in the present paper to deduce the 
position of these particles. It is the purpose of this paper to compare the time 
evolution of the ring current position with the observed temperature and electron 
density distribution in the topside ionosphere. An identificable SAR-arc repre- 
sents only that portion of emission of atomic oxygen caused by excitation of the 
electrons in the high energy tail of the thermal distribution, which is discernible 
above the background of emission due to ionosphere recombination. The SAR- 
arc therefore is a limited manifestation of the more general phenomenon of the 
region of elevated electron temperatures and depressed electron densities (see 
Figure 3 of this paper and Figures 8 and 9 of Norton and Findlay 1969, Roble 
et al. , 1971). The theory of Cole (1965) relied upon SAR-arc intensity measure- 
ments only for the purpose of estimating the electron temperature T^ . Since 
then, this parameter has become available from direct measurements and it is 
evident that this theory applies, a priori, equally well to the broad region of 
elevated T^ and depressed Ug . The exponential dependence of X 6300 A emission 
together with its linear dependence on n^ and the atomic oxygen number density 
n(0) must explain the variability of X 6300 A emission as a function of latitude. 
The most dominant of these parameters in this consideration is Tg , 

It is clear that topside T^ and Ug data distinguish two substantially different 
regions. One is associated with auroral bombardment and field aligned electric 
currents where there are correlated increases of both Tg and ng (see Findlay et 
al. , 1969 and Figure 3 of this paper). The other region equatorwards of the 
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auroral region is one of increased Tg and generally depressed tig (see Norton 
and Findlay 1969, Roble et al. , 1971, Figure 3 of this paper). This is the region 
over which there is heat conduction from the magnetosphere i It is with the evo- 
lution in time during magnetic siofms of'this latter region that the present paper 
is concerned. 

Cole (1965) suggested heat conduction to the ionosphere as a major sink of 
energy of energetic trapped particles. He proposed that as the trapped particles 
enter the inner magnetosphere coulomb collisions occur with the thermal plasma 
leading to a high background plasma temperature of approximately lev. Particle- 
wave interactions have also been su^ested as a mechanism to couple the energy 
of the ring current to the thermal plasma of the magnetosphere {Cornwall et al. , 
1971). This mechanism has recently been treated in a general fashion (Lin and 
Parks 1974). Conduction of heat occurs in the thermal plasma along the geo- 
magnetic field to the lower ionosphere and neutral atmosphere and causes an 
elevation of the electron temperature there. Recently Newton et al. , (1974) 
have suggested that enhanced electron temperatures in these regions may cause 
vibrational excitation of Nj which in turn can enhance the recombination rate in 
the ionosphere and cause some decrease of electron densities there. If the 
electron temperature is sufficiently great, emission at X 6300 A from atomic 
oxygen may appear above the airglow background. However, it has been shown 
that excitation of atomic oxygen emissions accounts for only a very small fraction 
of the total energy loss of the electrons, most of which is lost through other 
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processes (Cole 1965) determined by the electron temperature. Therefore the 
distribution of electron temperature in the upper thermosphere provides the best 
means of examining the heat flow from the magnetosphere to the ionosphere dur- 
ing storms. 

At the time of writing of his original paper Cole (1965) had only airglow 
measurements at his disposal for this purpose. Satellite measurements of 
electron density and temperature (Norton and Findlay 1969) and of ion densities 
and temperatures (Chandra et al. , 1971) on field lines above the SAR-arc have 
been reported for the storm of 29 September 1967. These data are consistent 
with the conduction theory (see also Roble et al. , 1971). Findlay (1969) reported 
that a narrow region of high electron temperature and a corresponding trough 
in electron density was a consistent feature found at altitudes of about 1200 km 
above observed SAR-arcs and that this feature moved in unison with the SAR- 
arcs through the storm period 31 October - 2 November 1968. 

Using these data Roble et al. , (1971) calculated 6300 A emissions on the 
basis of heat conduction. They then compared the results to observations and 
foxmd agreement, hi this paper we use the electron temperature distribution 
at 1200 km as identification of the position of the SAR-arcs and the aforementioned 
geomagnetic fluctuation data as an indication of the position of the ring current. 
The evolution in U. T. of electron temperature distributions and the spatial dis - 
tribution of geomagnetic fluctuations are studied here. These data in combination 
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support the theory of formation of the ring current (Cole 1964) and dissipation 
of its energy in the ionosphere (Cole 1965; Findlay and Cole 1970). 

OBSERVATIONS 

The storm of October 28 - November 2, 1968 is selected for detailed study 
because of the amount of data available in the form of satellite observations of 
n^ and T^ and ground observation of SAR-arcs at midlatitudes and of surface 
magnetic distrubance data. 

Figure 2 shows the variation of Kp through the period October 28 - November 
2, 1968, the times of satellite passes and the intervals of observations of SAR- 
arcs. Figure 2 also shows how the position of the inner edge of the relatively 
disturbed magnetosphere varies during the storm period as determined from the 
latitude profiles of magnetic fluctuations. This position is determined from 
plots of magnetic fluctuations as a function of invariant latitude. In the diagrams 
in this paper we have chosen invariant latitude rather than L~value as coordin- 
ate because the satellite data is more uniformly distributed in the former co- 
ordinate. The magnetic data were derived from the 2. 5 minute digitized ground 
based magnetic observations in the manner referenced above. The position of 
the inner edge of the magnetospheric proton fluxes is determined by drawing 
two straight lines (see Figures 4-9) one through the midlatitude disturbance 
and one along the steep rise of disturbance at higher latitudes. This rough 
procedure determines a point of intersection which, to within an accuracy of a 
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few degrees, iadicates the latitude at which high-latitude disturbance 
commences. 

We associate this surface latitude with the inner edge of observed proton 
fluxes in the magnetosphere on the same L-shell (Cole 1970 a,b, see Figure 1, 
this paper). Figure 3 shows the latitudinal distribution of electron temperature 
and density observed for the'satellite passes indicated on Figure 2. There are 
ground based observations of the SAR-arc for each pass except numbers 640 
and 652. For those two passes (as well as for the other) Roble et al. , (1971) 
have calculated the position of the SAR-arc. The calculations and observation, 
agree for the passes 662, 674 and 687. The position of the SAR-arc as indi- 
cated by the electron temperature distribution (hatched areas in Figure 2) in 
passes numbers 640 and 652 is poleward of ground observatories, thus photo- 
metric confirmation of the position of the SAR-arc was not possible. Figures 
4-8 combine plots of , Tg, and magnetic fluctuation profiles for five satellite 
passes approximately 24 hours apart. The magnetic disturbance profiles were 
determined in each case for the five hour period prior to the satellite pass and 
averaged over all longitudes. 

It can be seen that for passes 640, 652 and 662 a good correlation exists 
between the , T^ features and the magnetic profile which as noted earlier 
(see Figure 1), is an indication of the position of the ring current. However, 
such a good correlation is not found for passes 674 and 687. 
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In Figures 9 and 10 the satellite data for passes 674 and 687 are replotted 
against magnetic profiles taken during the five hour period in each of the major 
bursts in activity as indicated by Kp prior to the satellite pass concerned. For 
pass 674 we use the magnetic profile for the period 1530 UT ± 2, 5 hours on 
October 31. For pass 687 we use the magnetic profile for the period 1530 UT 
± 2. 5 hours on November 1. It is seen that good correlation is now also estab- 
lished for passes 674 and 687. In the case of pass 687 good correlation has also 
been found using several magnetic profiles taken during the active period from 
1200 hours UT on November 1 to 0600 hours UT on November 2; whereas, cor- 
relation between Ug, Tg features from pass 687 and magnetic profiles takeh prior 
to the increase in Kp at 0800 UT on November 1 was found to be poor. Later 
on in November 5 and 6, the daily average of Ap subsided to very low values 
6 and 7) and during this time the n^, Tg profiles returned to normal, but on 
November 9 there was a small burst of magnetic activity. This burst had the 
effect (on pass 746) qualitatively similar to that described on passes 650, 652 
and 662 and fitted the physical model tested in this paper. These observations 
can be summarized by saying the the n^ , Tg features correlate both in position 
and time with a previous burst of new magnetic activity. The persistence of the 
Tg profile is then explained by the lifetime of the particles injected deep into the 
magnetosphere at the time of the burst of magnetic activity (Cole 1964, 1965). 

Other general features in Ug, Tg distributions are evident in the data. For 
example note in pass 687 the broad region of high Tg and depressed ng extending 
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from ~45° to ~ 60°. In addition, there are observed anticorrelated spatial 
changes in Tg and in n^ in this region. This is distinct from the region in the 
vicinity of A 70° where both and Tg increase and which has previously been 
identified with auroral precipitation (Findlay et al. , 1969, Norton & Findlay 
1969). We identify the former region as the region of heat conduction from the 
magnetosphere which contains a similarly broad distribution of particles (for 
example see Figure 2 of Frank 1967 and Figure 4 of Cole 1970 b). 

DISCUSSION 

The foregoing observations indicate that for all conditions of magnetic dis- 
turbance the location in L-space of heat conduction from the ring current may be 
determined from the position of the region of elevated Tg and depressed Ug . It 
is important to reiterate (Cole 1965, 1970 a. Cole et al. , 1972) that the sink of 
energy for the ring current in the ionosphere is characterized by elevated Tg and 
that the SAR-arc as observed from the ground is only that portion of the region 
in which the value of Tg is sufficiently high that emission of 6300 A photons re- 
sulting from thermal exitation is detectable above the general background inten- 
sity in this line. This paper shows that the region of heat flow from the ring 
current to the ionosphere often covers a region two to three times as broad as 
any SAR-arc as detected from the ground. Thus, for example in Figure 3, the 
width of the region of elevated Tg is about 1500 km (from L *=» 2 to « 4) whereas 
the SAR-arc as observed from the groimd is only about 500km wide (from 
L 2. 42 to 2. 76). Evidence hearing on this point can be found in earlier papers 
(Norton and Findlay 1969, Roble et al. , 1971). 
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Theories of degradation of energy of the ring current must ejqjlain this 
broad region of elevated Tg. The theory of Cornwall et al. , (1971) is restricted 
in application to a narrow region (or order 500 km) just equatorward of the L* 
value of the plasmapause. 

The theory of Cole (1965) is not restricted in this way and predicts flow of 
heat into the ionosphere from the whole body of the ring current. The existence 
of cold plasma (a few ev) on field lines connected to the ring current is clear 
from the observations (Figure 3). The theory of Cornwall et al. , (1971) re- 
quires that the density in the equatorial plane be of the order of 100 cm"^ before 
their mechanism is effective in heating the electron gas. Coulomb collisions, 
however, as employed in the theory of Cole (1965) continue beyond the region of 
equatorial plasmapause where the densities of cold plasma may be less. Ener- 
getic particles with pitch angles at the equatorial plane less than about 45° will 
exchange energy by Coulomb collisions mostly near their mirror points where 
the background electron density may be up to 3 x 10^ cm'^ (Cole 1974, to be 
published). These particles will provide heat flow to the ionosphere at L-values 
beyond the plasmapause. 

We have demonstrated in the foregoing the spatial relationship of magnetic 
fluctuations and Tg profiles. Note the multiple peaks in Tg which may be related 
to the multiple bursts of magnetic activity. At this stage it should be mentitmed 
that multiple SAR-arcs have been reported (Roach and Roach 1963). An ex- 
amination of a wide range of data shows that in general the "pure" conduction 
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region is well separated from the precipitation region (Norton and Findlay 
1969). These properties of the profiles persist even at quiet times and sup- 
ports the contention that SAR-arcs exist at quiet times at high L-values (Cole 
1965). Schaeffer and Jacka (1971) have reported observations of X6200A 
skyglow that may be interpreted in this manner. In this connection it is 
important to note that the data on Tg and ng indicate no difference between 
the conduction region and the more limited region of the SAR-arc as observed 
from the groimd. 

Finally, within the resolution of the data available it is apparent that the 
conduction region featuring Ug and Tg appears immediately following each burst 
of magnetic activity. Moreover the data shows that the conduction region per- 
sists through the lifetime of the ring current. For example data on November 
14, 1968 at 1200 UT when Dst is approaching zero the conduction region is still 
discernible. 
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Figure 1, Schematic representation of the relation between ring current particles 
plasmapause and prior magnetic activity. (After Cole 1970 a,b) 
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Figure 2. Histogram of variation of Kp during the period October 28 ^ November 2, 1968 
together with variation of position of inner edge of relatively disturbed magnetosphere 
(dashed line) during this period. 
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Figure 3. Composite of all electron density and temperature pro- 
files versus magnetic invariant latitude for passes analyzed in the 
period of magnetic activity. 
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Plots of magnetic disturbance as function of invariant latitude together with topside and 
The magnetic disturbance relates to the average over the previous five hours. 
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Figure 5, Plots of magnetic disturbance as function of invariant latitude together with topside and n^* 
The magnetic disturbance relates to the average over the previous five hours. 
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Plots of magnetic disturbance as function of invariant latitude together with topside Tg and Ug, 
The magnetic disturbance relates to the average over the previous five hours. 
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Figure 9. Plots of topside Tg and Ug against invariant latitude for pass 674 at 8 hrs. UT November 1 together 
with magnetic disturbance for prior burst of activity during period 1530 ±2.5 UT on Oct. 31. 
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Figure 10. Plots of topside and against invariant latitude for pass 687 at 1000 hrs. UT on November 2 
together with magnetic disturbance for prior burst of activity during period 1530 ±2. 5 UT on November 1. 
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